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Figure 1: Illustration of the quantum spin liquid ground state. 

1 Introduction 

In these lectures we shall derive and discuss a set of exact eigenstates of the 
Haldane-Shastry jj], ^ model, a realization of the spin-1/2 Heisenberg chain in 
which the quantum-disordered spin liquid ground state and the neutral, spin-1/2 
excitations of such systems are particularly easy to understand. This behavior 
is not unique to the model, and in particular occurs in the Bethe solution of 
the near-neighbor Heisenberg chain, where it was first discovered || ||, but it is 
more accessible in this form. The model also makes the relationship of the spin 
chain to the fractional quantum Hall effect transparent ^ . The key results are 
these: 

1 . The ground state has the same functional form as the fractional quantum 
Hall ground state. 

2. This ground state has quantum disorder and has the same relation to 
the antiferromagnetically ordered state that a quantum liquid has to a 
conventional crystal. 

3. The elementary excitations of this state are spin-1/2 particles, spinons, 
and not spin waves, which are the elementary excitations of an ordered 
antiferromagnet. Spinons have 

(a) A relativistic band structure with a Dirac point at momentum ±7r/2. 

(b) 1/2 fractional statistics. 

(c) N/2 allowed momenta rather than N. 

4. The ground state is effectively 2-fold degenerate. The missing N/2 states 
of the spinon are excitations of the other ground state. 

5. The Hamiltonian has a factorization which may be construed as super- 
symmetric. 
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2 Haldane-Shastry Hamiltonian 



Let a lattice of N sites be wrapped onto a unit circle, as shown in Fig. 2, so 
that each site may be expressed as a complex number z a satisfying 



1 = 



(1) 



and let each site possess a single unpaired half-integral spin with its correspond- 
ing spin operator S a . The Haldane-Shastry Hamiltonian is then given by 



This Hamiltonian is translationally invariant, satisfies 

N 

[Hhs,S} = S = Y,So 

a 

and possesses the special internal symmetry 



IHhs,A] = Q a = ^ 



(S a X Sp) 



(2) 



(3) 



(4) 



Proof 

We proceed by applying the commutation relations 

[S 1 ,(S 1 -S 2 )]=i{S 1 xS 2 ) 



(5) 



the vector identity 



[(S 1 xS 2 ),(S 1 -S 2 )] = -(S 1 -S 2 ) 



a x (b x c) = (a ■ c)b - (a -b) c 



the fact that 

Zj + Zk Z£ + Zk 



z, - Zk zt — Zk 



ZjZkZe 



I Zj - ze 2 (zj - z k )(zi - Zk){zj ~ zt) 
which is totally antisymmetric under permutation of j, k, and I, and 

Zj + Zk 1 



j=£k z i Zk ' z i Zk ' 



\ z a (z a + 1) _ „ 



(6) 



(7) 



(8) 



(9) 
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Figure 2: Left: Illustration of Haldane-Shastry model. Right: Magnetic corre- 
lation function defined by Eq. (28). 



We have 



EE 

j^k a^fi 



Zj + Zk 



Zj - z k \z a - zp\ 



-[{S 3 xS k )AS a -Sp)} 



Zj + Zk 



ittt z j- z >°\ z j- z t\ 



+ i 



zj + z k 



Zj Zk I Zj Zk 



(Sj ■ Sk) Si - (Si ■ S k ) Sj 



^(Sj-S k ) = .□ 



(10) 



It follows from the commutation relations 

[S a , S b ] = i e abc S c [S a , A b ] = i e abc A c , 

that Hhs> S 2 , and (A • S) all commute with each other. 



(11) 



3 Ground State 

Let us imagine the spin system to be a 1-dimensional string of boxes populated 
by hard-core bosons, the J. spin state corresponding to an empty box and the j 
spin state corresponding to an occupied one. The ground state wavefunction is 
a rule by which a complex number is assigned to each boson configuration. The 
total number of bosons is conserved, as it is physically the same thing as the 
eigenvalue of S z . Let N be even and let zi, denote the locations of N/2 
bosons, the number appropriate for a spin singlet. Then the Haldane-Shastry 
ground state is 
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N/2 



N/2 



...,z N/2 ) = Y[(zj - Yi 



(12) 



j<k 



3=1 



Its energy is 



H U , *>=- (|y) ( V • -| ) : 



(13) 



Proof 



We begin by observing that [S+Sp &](zi, z N / 2 ) is identically zero unless one 
of the arguments z\,...,z N / 2 equals z a . We have 



N/2 N 



AT/2 JV-2 r JV 

ee{e 

j=l £=0 ^f3^j 



N/2 

E 

3 = 1 



(N-l)(N-5) N-3 2 d 

— Zi Za 



12 



2 3 dzi 



(N(N- l)(iV-5) 2z. 



24 



7— 



l 

+ 2 



AT/2 

E 



N/2 

(,-,)(,-^ + S ra 



2z? 



|*(^i,...,^ /2 ) 



G 



f N(N -1){N -5) N(N — 3) N 1 
\ 24 4 ~ -* + 2 



AT/2 



|*(2;i,...,^ /2 ) 



Here we have used the fact that 

1 



(14) 



(15) 



|Z Q -Z/3| 2 (z a - Zfj) 2 

i.e. that the Hamiltonian is effectively analytic in the spin coordinates, and the 
sums (with sites labeled so that zn = 1) 



N-i 1 

^ Ua-l| 2 = 



A^ 2 - 1 



AT-l 2 



E 



r! (*<, - i) 

N-l 

E 



JV-l 9 

z 



12 

(N - l)(N — 5) 
12 

AT-3 



r! (^ - i) 2 

(V-l 
a=l 

^^ -i)=...=^^-ir 



JV-l 



(16) 
(17) 
(18) 
(19) 
(20) 



a=l 



a=l 



+ 



(z a - zp)(z a - z 7 ) (zp - z a )(zp - Zj) (Zj - Z a )(Zj - Zp) 
which arc worked out in Appendix A. We also have 



= 1 



, N QZ QZ \ 

I \ - J a J )3 I 



. . .,z N/2 ) 



(21) 
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f N(N 2 - 1) ^ 1 It/ 



(22) 



This completes the proof, since 



«»-^)'{Si^ + Si^}- D ,23) 



Important Properties 
1. Reality: Since z a lies on the unit circle we have 

N/2 N/2 

**( Zl ,...,z N/2 ) = l[(z*-zt) 2 l[ 

j<k j 



N/2 N/2 

= tl{z k -z j ) 2 l[z]- N = *{z 1 ,...,z N/2 ) . (24) 

3<k j 

Thus * is real despite being a polynomial in the complex variables Zj . 

2. Translational Invariance: ^ is translated one lattice spacing by multi- 
plying each of its arguments by z = exp(i27r/iV). Since it is a homogeneous 
polynomial of degree (N/2) 2 we have 

*(ziz,...,z N/2 z) = exp(i/V7r/2) *(z lt . .. ,z N/2 ) . (25) 

The crystal momentum of the state, i.e. the phase it acquires under trans- 
lation, is thus or 7r, depending on the value of the even integer N. 

3. Spin Rotational Invariance: To prove is a spin singlet it suffices 
to show that it is an eigenstate of S z with eigenvalue zero and that it is 
destroyed by the spin lowering operator S~ . The former is true for any 
wavefunction in which the number of Zj is constrained to N/2. For the 
latter we have 

N 

[S~^]{z 2 ,...,z N/2 ) = ^2^(z a ,z 2 ,...,z N / 2 ) 

n=1 



a=l 

N-l ( N , e 
t=l ^a=l ' 1 



Z!,z 2 ,...,z N/2 ) =0 , (26) 



8 



since 



N 

Y,zi = N5 m (mod N) . (27) 

a=l 

This implies that the wavefunction is the same with the roles of | and 
J. reversed or, more generally, with the quantization axis taken to be an 
arbitrary direction in spin space. 

4. Quantum Disorder: In Fig. 2 we plot the spin-spin correlation function 

<S a ■ Sp > 



= 3N E Z3 ,..., ZN/2 \^^ a ,z (3 ,z 3 ,... 1 z N/2 )\ 2 3 

2 l 2 ' E Z2 „,„z N/2 \V(z a ,z 2 ,...,z N/2 )\ 2 4 • [ > 

The convergence of this function to zero as \z a — zp\ ~ + oo shows that |\P > 
has no long-range order and is a spin liquid. The fall-off is slow, however, 
and this is important. Strongly-disordered spin liquids, i.e. with exponen- 
tial decay of correlations on a length scale £, are easy to construct when 
the spin in the unit cell is integral. They have an energy gap A = hv/£, 
where v is the spin- wave velocity of a nearby ordered state. But strongly- 
disordered spin liquids cannot be stabilized with short-range interactions 
when the spin in the unit cell is half-integral. The excitation spectrum in 
this case is always gapless || . 

5. Factorizability: The fact that > is a product of pair factors makes 
a number of its properties easy to calculate by semi-classical monte-carlo 
techniques. Let us, for example, consider Eq. (EST) . Writing 



|*(zi, z N/2 )\ 2 = exp 



[zi, -,z N / 2 ) 



(f>{zi ) ...,z N / 2 ) = -4y^lnlzj -Zfc| , (29) 
j<k 

we see that the summand is the Boltzmann factor of an equivalent finite- 
temperature classical lattice gas, and that we are computing the joint 
probability for two of these particles to reside at sites z a and zp. This 
may be done by generating a time sequence of configurations using a rule 
that obeys detailed balance, and then simply counting how many times the 
sites z a and zp are simultaneously occupied. The simplest such algorithm 
is the following. Let the current configuration be Z\, 2jv/2- 



(a) Loop on particles j. 

(b) For this particle, roll the dice to choose a direction. Compute 

z'j = Zj exp(±i2n/N) 

depending on the outcome. 

(c) Compute 



(30) 



/=n 



Zj Zk 



(31) 



(d) If / > 1 update Zj to z'j. 

(e) If / < 1, roll the dice to generate a real number x between and 1. 
Update zj to z'j if x < f but do nothing otherwise. 

Fig. 2 was generated using this algorithm. 

6. Degeneracy: The Haldane-Shastry ground state is not degenerate, but 
it is nearly so. The alternate ground state is 



N/2 



N/2 



y , (z 1 ,...,z N/2 )=Y[(z j -z k f i~n 

j<k L j=i 

It has crystal momentum tt greater than that of ^ and has energy 

J\ / 7 



(32) 



Hhs\^>=-J 



v -- ) |*'> 



(33) 



It is equivalent to the original vacuum plus a pair of spinons excited out 
of the vacuum into a singlet with total momentum n. 



4 Spinons 



Let the number of sites TV be odd and let 

M M 



M 



*a(Zl, -,Zm) = Y[(z a - Zj) "[[(Zj - Z k ) 2 Y[ 
3 j<k j 



(34) 



where M = (N — l)/2. This is a j spin on site a surrounded by an otherwise 
featureless singlet sea. We have 
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Figure 3: Spin and charge profiles of the localized spinon \^> a > defined by Eq. 
(34). The dotted lines are a guide to the eye. 



JV 



per Eq. (Eq). The combination of these states given by 



N 



Zl, 



, zm) , 



(35) 



(36) 



with < m < (N — l)/2, is an exact eigenstate of the Hamiltonian with eigen- 
value 



n H s\^m>={-J(^)(N--) + -( — ) 2 m{—^- -m)}\* m > . (37) 



Proof 



Following Haldane Q we consider a wavefunction of the general form 



M 



M 



V{zi,...,zm) = ®{zx,.-.,z M ) Y[(zj -z k ) 2 Y[; 

j<k j 



(38) 



where M = (N — l)/2 and $ is a homogeneous symmetric polynomial of degree 
less than N — 2M + 2. This latter condition causes f to be a polynomial of 
degree less than N + 1 in each of its variables zj, and thus allows the Taylor 
expansion technique used for the ground state to be applied. Doing so, we find 
that 
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provided that <& satisfies 



3 j^k Zj Zk dZj 

Let us now consider the polynomial 



E 



9$ 



(39) 



(40) 



M 



M 



$A(zi, . . . ,2lf) = \\{ Z A - Zj) = ^ z™P m (zi, . . .,Z M ) 
j m— 

where za is not necessarily a lattice site. We have 



(41) 



m 



2\^** c)z- 

3 



2^ r, _ z, flz, [ 9 ^ 



J ^fc ^ - Zk dz i 



= *Y 



z 2 



+ 



(Zj - ZA)(Zj - Zk) {Zk ~ ZA){Zk - Zj) 



M 



N - 3 z j 



Zj - Z A 



= { 



{ M(M-l)-z Ad ^ 



d 2 N-3 



M-z A 



3za 



<$>A 



M 



m=0 



(42) 



The proof is completed by multiplying both sides of this equation by (z* A ) n and 
then summing on lattice sites za- □ 

The state | $ m > is a propagating { spinon with crystal momentum 



q=-N-—(m+-) (mod 2tt) , 



(43) 



per the definition 
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-7T -7T/2 7T/2 7T -TT -Tt/2 7r/2 TT 



Figure 4: Left: Spinon dispersion given by Eq. (46). Right: Allowed values of 
q for adjacent odd N. 



* m (ziz, . . . , z M z) = exp(iq) V m (zi,..., z M ) ■ 
Rewriting the eigenvalue as 



we obtain the dispersion relation 



(mod 7r) 



(44) 



(45) 



(46) 



plotted in Fig. 4. Note that the momenta available to the spinon span only the 
inner or outer half of the Brillouin zone, depending on whether N — 1 is divisible 
by 4. The spinon dispersion at low energies is linear in q with a velocity 



^spinon — ~^ ~jT~ 



(47) 



where b — 2n/N is the bond length. 

The existence of spinons is an automatic consequence of quantum disorder 
whenever the total spin in the unit cell is half-integral. The spin liquid ground 
state must be a singlet because otherwise it is ferromagnetic. A singlet is possible 
only if the total number of lattice sites N is even. If N is odd then the total 
spin can be no less than 1/2. But since there can be no physical difference 
between even and odd in the limit of large N, the system must have had a 
neutral spin- 1/2 excitation to begin with. 

The ground state of the odd-N spin chain is 4-fold degenerate and is given 
by \^ m > f° r m — and (N — l)/2 and their f counterparts. This corresponds 
physically to a "left-over" spinon with momentum ±7r. 

Spinons maintain their identity when more than one of them is present. 
When N is even, for example, the states 
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N N 



a=l 0=1 



where 



M MM 
* Q/3 (zi, ...,z M ) = Y[(z a - Zj){z - Zj) Y[ ( z j - Zkf J| Z 3 ' 

3 3<k j 



with M = N/2 - 1 and 



N/2 



frnn(z) = ^ dp/ - ^-a N/2 ; 



N/2 



(48) 



(49) 



m — n + 2£ \ - 
a ^2^ + m-n-l/2]^ o afe (a0 = 1} 

with to > n, are eigenstates of the Hamiltonian with eigenvalue 

„, ,t f w^x/,, 19 24 N J,2tt n2 

W **l*™ >= - J ( 24 )(7V -N + N^ )+ 2 ( iV } 



(50) 



TO (__l_ m ) +n (__l_ n ) 



to — n 




I 


2 







|*ron> • (51) 



Proof 

Following the procedure we use for one spinon we take $ to be a superposition 
of states of the form 

M 



®AB = ]J(ZA Zj)(z B - Zj 



M M 

= Yl Z A Z B Pm(zi, ...,Z M ) Pn(zi, ...,Z M ) 

m— n— 

where za and z# are not necessarily lattice sites. We find that 



(52) 



2 \Y j dz) f£zj-z k dzj J 2 ^ 



dzj 
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-(— ) 

V 2 ' 



1" Z B^. 

dz A oz B 



2M 2 - M(N - 2) 



>$AB 



m) + n (~^ — 1 — ti) 



m-n Z A + ZB t.m.np p 

2 2^ - Z B 



and thus that 



j2 4z j dq> ™ 

Zj Zf~ dzj 



M 



iV-3f d>S> mn 
2 4^ 92, 



N _/V 771 — 77- 

m (y - 1 - m) + n( y - 1 - n) H — ^$ TO n 



e=o 



for m > n, where 



N N 

= Z)( 2; a) m (^) n *a/3 = ^« • 
a=l /3=1 

In obtaining this last expression we have used the identity 

m — n 

' + y .(x m y n - x n y m ) = 2^2 x m ~ e y n+e - (x m y n + x n y m ) 



x-y 



The solution of Eq. (|4C|) is then 

n 



N N m — n 
A = m {~^ — 1 — m ) + n \~~2 — 1 ~ n ) ^ — 
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\og(z) 





7T —TT 





\og(z) 



Figure 5: Left: Real part of spinon pair internal wave function f mn (z) denned 
by Eq. (50) for the case of N = 100 and m - n = 2,4, 10 (solid) and 50 
(dots). Right: The "adjoint" function f mn (z) defined in Problem 7. 



where the coefficients ag are given by Eq. ( p0| ) . Such a simple solution is possible 
because the matrix to which Eq. (|54|) corresponds is lower triangular, i.e. takes 
the form 



Matrix = 



E 

Uio E x 

W20 «21 E 2 

«30 «31 «32 £3 



(58) 



The eigenvalues of such a matrix are its diagonal elements, and the correspond- 
ing eigenvectors are generated by recursion. It should be noted that the upper 
bound on the sum in Eq. ( |50| ) is flexible, as $ mrl is identically zero unless 
< m, n < M (mod N). We have chosen the largest possible value so as to 
optimize the smoothness and short-rangeness of f m n{z)- n 



It may be seen in Fig. 5 that f m n(z a / 'zp) exhibits a scattering resonance, 
an enhancement when z* a z$ ~ 1, indicating that the spinons attract each other. 
This attractive force may also be inferred from the energy eigenvalue if we 
rewrite it as 



H HS \* m n>={-J(^)(N + ^)+E qi +E q2 +V qi - q2 }\* mn > , (59) 
where 

T 27T . 1. 7T 2ir , 1. . „. 

Qi = (m+-) q 2 = (n+-) (60) 

H 2 N 2' 2 N 2 y ' 

are the spinon momenta, E q is defined as in Eq. ([l6"|), and 
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Vo = -J- 



N 



(61) 



Note that this potential vanishes as N — > oo, as expected for particles that 
interact only when they are close together, and that the total crystal momentum, 
as defined by Eq. (E3) , is 



q= -{N -2) + qi + q 2 (mod27r) , 



(62) 



a value greater by it than that of the ground state when qi = q 2 = 0. 

For every triplet state |\& mn > defined as in Eq. ( ^8|) there is a corresponding 
singlet A z S + \^ mn >, where A is defined as in Eq. (||), with exactly the same 
energy eigenvalue. For example, for the case of m = N/2 — 1 and n = we have 



N/2 



N/2 



[h z S + ^ mn ](z 1 ,...,z N/2 ) = —N 2 Y[( Zj ~z k f 

j<k 1 j=l 

which may be seen to be the alternate ground state defined in Eq. (|3 



(63) 



Proof 



We begin by using the sum rules 

M N—M 



n z j n % = 1 



M N-M 
z 3 Y[( z 3 ~ z k) II ( Z 3 ~ ? ' fc ) = N ' 



k 



to rewrite the wavefunction in terms of the j spin locations r\j , per 



N/2-1 

*mn(zi, -,ZN/2-l) = f| ( 
3<k 



We then have 



[5 + * m „](zi, Z N / 2 ) = ^2 ^mn(Vl, ■■; VN/2, Z a ) 



N/2 



j<k 



N/2 



= N\{{r l] -r lk f 1 + n^ 2 = Nj{{ Zj -z k f 1 + 



3 = 1 



N/2 



j<k 



N/2 



3=1 



and thus 



(64) 



N/2-1 N/2+1 

n fe-^) 2 • ( 65 ) 

j<k 



(66) 



[A 2 5+*„ m ](z 1 ,...,z Ar/2 ) 
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^£(7^)[^S+viw](zi,...,^/2) 



2^ z a - Z/3 



| £ £( Zj •■■> z J-l>2/3 5 Zj+1 5 ZjV/ 2 ) 



2 •« — ' ' — ' ' zh — zr 

3 0& 3 ' 



N/2 N/2 N/2 



JV/2 r N/2 . 

i<fe L j=i ■ 



□ 



(67) 



The singlet has no simple exact representation in terms of \^ a /3 > but is rea- 
sonably approximated by 



A z S+\* mn > 




That it is an energy eigenstate follows from the conservation of A and the non- 
commutativity of A and S implicit in Eqs. (0). That it is a singlet follows 
from 

A z \V mn >= \^~y~ -m-nj|* ron > , (69) 

i.e. that Y&mn > is an eigenstate of A 2 , an important result we shall revisit. This 
implies that the spin-2 representation contained in the 9 states h^S v \ i ^ mn > 
{n,v = 1,2,3) must be identically zero and the spin-1 representation must be 

just \^ mn > itself. 

Spinons are semions, i.e. particles obeying 1/2 fractional statistics. Since 
the 2-spinon wavefunction 

^ab = Y[{zj - z A )(zj - z B ) \\_{zj - z k f , (70) 

3 3<k 

where za and zb are not necessarily lattice sites, has the property 
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^*ab( z i, •••) Zn/2-i) = {zazb) 1 n/2 ^ab(zi, ZjV/2-l) 



(71) 



the Berry phase vector potential for adiabatic motion of spinon A in the presence 
of B is 



i d . d , . , 

<1PAB\ZAS WAB> + <ZA7 WAB\WAB> 

OZ A OZ A 



I <i>AB\4>AB> 



1, A\ 



(72) 



The phase to "exchange" the spinons by moving A all the way around the loop 
is thus 



A<j) 



IT 

±-i 
2 



(mod 2ir) 



(73) 



1(1 _ N \ dz A 

2 1 2 ' z A 

This number is or 7r for bosons or fermions. Fractional statistics is actually the 
long-range force between the spinons manifested in the resonant enhancements 
of Fig. 5 and the potential V qi - q2 in Eq. (59), and has nothing to do with the 
symmetry or antisymmetry of \^ a /3 > under interchange of a and (3. It does 
have to do with state-counting. The number of states available to £ I spinons, 
determined by counting the number of distinct symmetric polynomials of the 
form 



{N-i)/2 



<5> ZAl ,..., ZAe (z 1 ,...,z (N _ e)/2 ) = Yl {zj - z Al ) X ... X (zj - z Ae ) , (74) 



is 



A/?eroi=( iV/2 + ^/2 ) 

This is just halfway between the numbers 

TVf rroi = ( N/ . 2 ) A/? oso = ( N/2 . + £ ) 



(75) 



(76) 



likewise calculated assuming that the number of states available for one particle 
is N/2. 
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5 Annihilation Operators 

The operators 



2 £7' \z a - z/j 



annihilate the Haldane-Shastry ground state, i.e. satisfy 

n a \^>=o , 

for all a. 



(77) 



(78) 



Proof 



We have as before that [S+ Sp ^](zi, z N / 2 ) is zero unless one of the arguments 
z\, z N /2 equals z a , but in this case we do not sum over z a . Instead we have 



N 



N 



[ S a S p^]( Z a,Z2,...,Z N/2 ) = Y ^ ^( Z P^ Z 2,---, z N/2 



fj^a 



Za - Zfj 



(j^a 



Za - Zfj 



ZgZ(j{zp - ZgY 1 d 1 f^(Zg, Z2,..., Z N / 2 ) 



= — + 2^——\*(z a ,Z2,...,z N/2 ) . 



(79) 



However (1/2 + S , „)| V E' > is also identically zero unless one of the arguments 
Zi,...,Zn/2 equals z a - We thus have 



JV 

E 

fj^a 



+ S Z a ){\ + S^]( Za ,Z 2 , Z N/2 ) 
Z a ~ Zfj Z A 



= zC Z 1 - ^(Za,Z2,...,Z N/2 ) ■ 



(80) 



Subtracting these from each other we find that 

N 



r JV 



Zq - Zfj 



SZSp-2(\ + S' a ){\+S%) 
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|*>= 



(81) 



for all a. However since \^>> is a spin singlet the irreducible representations of 
the rotation group present in this operator must destroy |*> separately. The 
scalar component is identically zero. The vector component is 



N 



0? a 



Z a - ZfJ 



[i(S a xSp) + S }\*>=O 



(82) 



The rank-2 tensor component is the product of the two and therefore also zero. 
Since > is also its own time- reverse it must be destroyed by the time-reverse 
of the vector operator, i.e. 



E t^tt^" x 30 + sp] = - E T^rr^* x ^) + ^] • ( 83 ) 

The difference of these is the trivial operator S a x S, and their sum is 2f2 Q . □ 
These operators satisfy 



-On 



(84) 



and are related to A by 

]Tfi Q = A . (85) 

a 

They are not symmetries of the Hamiltonian but supercharges, for we have 



z a - 



16 



(86) 



exactly. 



Proof 
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Since 



ti a = J2 —?2—[i0 a x fy) + §p] - hi(S a x S) + 3\ , (87) 
^ z a - zp 2 



P^a 

per the previous discussion, we have 

= £^n a + |s.A + |(iV-l)S 2 . (88) 
In evaluating this expression we have used 

[i(S a x S 7 ) + 5 7 ] f = [z(S 7 x S a ) + S 7 ] , (89) 

^[i(^x^ Q ) + 5]-0 Q = ^[iS ? -(5 Q xO ct ) + S ? -0 Q ] = ^5-A , (90) 

a a 

and 

^[i(SxS a ) + S]-[i(S aX S) + Si 



2 

which follows from 



3 -Y,lS 2 -^S■(S a xS)} = llN-l}S 2 , (91) 



[i(5 7 x S Q ) + 5 7 ] • [i(5 Q x 5^) + ty] 



= -(S 7 x S a ) • (S a x Sp) + iS 7 • (S a x Sp) + z(5 7 x S a ) ■ Sp + S 7 • Sp 



= (1 + |)(£ 7 • S^) - (S 7 • S a )(S Q • S>) + 2z5 7 • (5 Q x Sp) 

= ^{S 7 -Sp + iS^(S a xSp)} . (92) 



The 2-spin sum is 



^ Z*)(z a - Zp) (Z a -Zj)(z a -Zp) 7 



P^a ^ a 7/v « V> a#/3# 7 
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Er^v^ E 



^7 



- Z 



2/3 



E z 7 (z 7 + Zp) g - 1 { Z/3 + 2 - , q r- 



0/7 



2 8 ^ |z Q - z^ 



The 3-spin sum is 



E 



g 7 ■ (Sg X 5,3 ) 

(z* - z*)(z a - zp) 



y a 7 



; E 



(z a - z 1 )(z a - zp) 



S a ■ (5 7 x Sp) 



z a (z 7 — Z/3) 



S a ■ (5 7 X Sp) 



- y 

2 ^ 

a#/3#7 



z a z 7 Z a Zp 



S a ■ (5 7 X 5/3 ) 



a//3# 7 



(5 Q x 5 7 ) • Sp = A • 5 . 



Note that in the last step we have used the identity 



{S a X 5 7 ) • {S a + 5 7 ) = 

Putting these results together, we find that 



EEEtt^ 

—^——^[-i{S a x 5 7 ) + 5 7 ] • [*(£„ x 5^) 



a p=£a j^ia 
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2 ir 1 ^r 1 {z* a - z*){z a - zp) 



[S 1 ■ Sp + iS 1 ■ (S a x Sip)] 



2 1£ i* 



■ 3 S a • 5/3 



zp\ 



N(N 2 + 5) 
16 



□ 



(96) 



Since < -fi Q |$ > is non- negative for any wavefunction |$ >, this provides an 
explicit demonstration that |W> is the true ground state. The annihilation op- 
erators and their equivalence to TLhs when squared and summed were originally 
discovered by Shastry R] . They are modeled after a similar set of operators dis- 
covered for the 2-dimensional chiral spin liquid, although there was a minus-sign 
error in the original paper which caused the operators to be mistakenly reported 
as scalars ||. They are lattice versions of the Knizhnik-Zamolodchikov opera- 
tors known from studies of the Cologero-Suthcrland model, the 1-dimensional 
Bose gas with inverse-square repulsions [l^] . 

6 Spin Current 



The operator A is a scaled spin current. Its action on the propagating spinon 
of Eq. (|36l), for example, is 



■ N- 1 

A 2 |* m >=<! ; m}\^ m > , 



which is proportional to the spinon velocity 



dE q 
dq 



2tt J f N - 1 



N 



(97) 



(98) 



Its action on the 2-spinon state given by Eq. ( p9|) is similarly the sum of the 
two spinon velocities. 



Proof 
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We have, with M = (N - l)/2, 



M 



[A z V A ](zi, ...,z m ) = \^2^2(^-^-)^a(zi, ...,Zj-i,zp,z j+ i, ...,z M ) 

Zj Zp 



3 0¥=j 



\ E E \ [E(^)(^ - h{* A ^ zm),z ] 

= {—^-zAq^I*a(z 1 ,...,z m ) , 



and thus 



A z \^ m >=J2(z A ) m A z \^ A >= ^-ro}|* r 



> . □ 



(99) 



(100) 



A more traditional description of this current may be constructed by interpo- 
lating the spin operators into the interstices by means of the formula 



a(z) = 

The Hamiltonian is then 



z N/2 _ z -N/2 

2iV 



E 



z + zp 

z - Zfj 



S(3 ■ 



(101) 



1 f. da . , da^dz 



N 

E 



S a ■ Sp 3 



2ni J dz dz z N '—f, \z„ — zg\ 



(102) 



where the integral is performed over the unit circle. We also have spin density 
and spin current density operators 



p{z) = —ia(z) x a(z) 



->, , 1 ( r da da 



which satisfy the continuity equation 



lim 



\z a - zp\- 



,P\ =0 . 



The zero-momentum component of this conserved current density is 

2iri J z 



(103) 



(104) 



(105) 
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7 Supersymmetry 



We shall now consider the generalization of the Haldane-Shastry Hamiltonian 
~~ ^ 5Z( c l s c /3s + c]j s c QS ) + ~(n a + n ) - ^n a np ~^\p (106) 

S " J 

first studied by Kuramoto and Yokoyama where 

p = n (! - c iAi c -i c ^ > ( io? ) 

is the Gutzwiller projector, and site occupation and spin operators are 

n a = C^Caf + C^CqX = o( C c»T C °4 + C^C a f ) 



S a = Yi^ Cal ~ C "1 Cq T) S a = ^( c i T C «T " c al c al) • (108) 

Thus each site can be f, |, or unoccupied - but not doubly occupied - and the 
hole can tunnel to nearby sites by means of the same inverse-square matrix ele- 
ment characterizing the spin exchange. This equivalence of the energy scales for 
magnetism and charge transport causes the Hamiltonian to be supersymmetric 
in the sense of Eq. ( jS6| ) and also in the more traditional one of commuting with 
electron or hole injection per 

[H K Y,Pc as P]=0 . (109) 



Proof 

It suffices to show that 



where 



[n af 3,(c a] + c m )] |v»=o , (no) 

H a p = P< S a ■ Sp ~ j ^2(cl s Cp s + Cp s C as ) 
^ s 

- ^n a np + i(n a + np - 1) j P (111) 

for the 9 configurations \ip > on sites a and j3 allowed by the projector P. 
Denoting the state with no electron on either site by |0> we have 
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Case 1: 



l^i>=cL T 4 T |0> (112) 



(CaT+C/Jt)Wa J 8|^l>= | (c qT + C/3 T ) | >= J(cJ ? -C^)|0> 



H Q/3 (c aT +c (3T )|?A 1 >=H q/ 3(c| 3T -c^ t )|0>= -(c+ t -c^ t )|0> 
Case 2: 

iV'2>=(ci T 4 i -44 T )|o> (H3) 

(c aT +C / 3 T )W aj8 >= -7j(c«T +C/5 T )|V2>= - ^ (c^ + C^) |0 > 
Ha/3(Ca T + C /3T )IV'2>= Wa/j(cJ A + C^)|0>= + C^)|0> 

Case 3: 

iv>3>=(4 T 4 i+c t ei 4 T )|o> (H4) 

(c qT +C /3T )W a /3|V'3>= ^ (c qT + C/3 T ) | V>3 > = ^(^i " C al)l°> 

W a/3 (c aT +c / 3 T )|V3>=W a/3 (4 i -c^)|0>= i(4 x -c^)|0> 
Case 4: 

|V4>=4 T cJ t |0> (115) 



(c« T + C^Wc^ |V>4>= ^(CaT + C/3t)IV'4> = 



H a (i{c a ] +C /3T )|V'4>= 
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• Cases 5 and 6: 

|^5>=c qT |0> (116) 

H q/3 (c qT + C^)\lf) 5 > = H a p\0>= ~|0> 

• Cases 7 and 8: 

|^ 7 >=c Qi |0> (117) 
H q/3 (c qT +c /3T )|V'7>= 

• Case 9: 

|V^9>=|0> (118) 

(CaT + C^Wc^l^g >= ~^( C «T + c /9t)I°>= 
HotpiCotf + C/3-[)\lp9>= . □ 

8 Holons 

Holons are charged, spin-0 elementary excitations of the Kuramoto-Yokoyama 
Hamiltonian made by removing an electron from the center of a spinon. The 
localized holon CaJ^Q >, where |\I> Q > is defined as in Eq. fl34|), is shown in Fig. 
6. It is a natural complement to the spinon from which it was made, for the two 
states differ only in the disposition of the central site, which does not fluctuate 
and is not involved in any way in the quantum number fractionalization. This 
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Figure 6: Spin and charge profiles of the localized holon c a J^/ a >. The dotted 
lines are a guide to the eye. 



state is an exact spin singlet by virtue of Eq. (|35|). Let N be odd and let 
M = (N — l)/2 as in Eq. (|36|). Then the propagating holon wavefunction 

M M M 

*% l ™(z 1 ,...,z M \h) = (h*) m l[(h-z j )l[(z j -z k fl[z j , (119) 

j j<k j 

where z\, Zyi denote the positions of the \ sites and h denotes the position 
of the empty site, all others being J., satisfies 

n K Y\^: lon > 



Is J ' 2?!\ 



N + 1 



\ V 24 A N 2 N V 2 
provided that -(N + l)/2 < m < 0. 



to) ^1* 



.holon 



> 



(120) 



Proof 



To simplify the notation let us negate the power of h in Eq. ( |119| ) so that the 
holon wavefunction becomes 



M 



M 



* n (g 1 ,...,z M \h) = h n Y[{h-z j )'[[( 



Zj - z k 



M 

n 



(121) 



j<k 



with < n < (N + l)/2 and M = (N - l)/2. The effect of the spin-exchange 
part of the Hamiltonian is the same as for the spinon. Taylor expanding as 
usual we obtain 
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N s + s~ 

a 13 * n ](z u ...,z M \h) 



\z a - zp\ 



1-N 2 



M 1 
E U. _ . 



24 ^ k \Zi-z k \ 



2 dh dh 2 J \ /i™ J v ; 

The charge- exchange terms also behave similarly. The operator c a ^cj^ gives 
zero unless the holon resides at z a . For this case we have 



N 

E 



z a - z \ 2 



N 1 

*„](zi, Z M \z a ) = ^ I Z |2*™( Zl ' z M\zp) 

\ z a Zfi\ 



E' y^ 1 zp(zf) - z a ) e d t \ ^„{zi,...,z M \z a ) 
f-^ £\\zB-z a \ 2 { dzj 

j3^a £=0 1 P 1 



Z" 



N 2 - 1 n{n - N) 



12 



N — 1 



z n+l d + l _ z n+2 9 2 \ f ^ n (zi,-,Z M \z a ) 



a dz a 2 



dzl 



(123) 



Summing on a we then obtain 



N Pr ir f P 

|-2q - Z/3| 



Et 



*„](zi, ...,z M \h) 



N 2 - 1 n(n - AT) 



12 



N-l 



-.n+l JZ. i 1 »."+2 Q2 \ f ^n(zi,...,Z M \h) 

dh + 2 dh 2 \\ h" 



(124) 



For the other charge-exchange channel we use the fact that \^ n > is the same 
written in the J, coordinates 771, ...,t]m by virtue of being a singlet. This gives 



n p t p 



*™](»7i, -,VM\h) 



N 2 - 1 n(n - TV) 



12 



+ 



iV-1 



— n 



8 1 <9 2 
h n+1 — + -h n+2 — 
dh + 4 



^n(?7l, ...,7?M|/t) 



(125) 
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It remains only to rewrite this expression in terms of the f coordinates zi, zm- 
We have 



h n+1 - 



d \ ^„{r]i, ...,r) M \h) \ ft 



9ft ft™ 



f = Vn(m,-,VM\h) 

J , ft - Vj 



= {^V-E^}*«(^'-^Mift) 

={(^-ft-l-}{^i*} , (126 ) 



and 



■ n+2 g2 f ^n(?7l,-,T?M|ft) ] = ^ 

9ft 2 1 ft" J 4^ 



ft 2 



(ft - 7]j)(h - rjk) 



*n(»7l, -,»?M|ft) 



U ft fJV-1 ^ ft ft 



5 ft - % I 2 ^ ft 



Z fe ft - T]j 



*n(»?l, -,»?M|ft) 



M 



N — 1 y-> ft 

ft" 



2 ^ h - z k 



+ 



(iV-l)(iV-5) 
12 



(TV - l)(AT-2) 



M 



a 2 l r*„(z 1 ,...,z M |ft) 



ft" 



which gives 



A r Pc aT c t P fn(n-l),„ J£, 

E[^^ F ^](^-.^ift) = {^-^ft" + E 



(127) 



ft"+ 2 
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— n h 



n+id_ 1 n+2 d 2 \ f ^ n (z 1 ,...,z M \h) 
dh 2 a/i 2 J 1 ft" 



For the Ising and electrostatic energies we have 



N 

D 

a =£f3 



CZ QZ 

ct {3 

\z a - zp\- 



:*n](zi, z M \h) 



(128) 



M ^ M 

E U, _ r ,|2 + E 



1 N(iV 2 - 1) 



and 



AT 
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V n (z 1 ,...,z M \h) (129) 



y 1 



1 13 

-(n a +n^) - -n a ^ - - 



$> n (zi, Z M \h) 



1 - TV 2 
24 



Adding these contributions together we obtain 

N 



E ii i2 p 5 « 5 /j + s « s l + E + + n p) 



\Z a - Z/3\ 



(130) 



- A n a n p - - 



N(N 2 -l) , iV+l x 
— +n(n -— ) 



P V n ( Zl ,...,z M \h) 

* n (zu...,z M \h) . □ (131) 



The propagating holon wavefunction differs from that of the spinon in being 
defined for all values of m, including those for which it is not an eigenstate of 
the TLky- For these other values of m we observe that 



N 



o< i<* rn |^p4 i p|*^ o ">| 2 <i , 



(132) 



i.e. that a supersymmetric rotation of the spinon, which is an eigenstate of 
H.KYi has a nonzero projection onto the corresponding holon. The physical 
meaning of this is that the holon can lose energy by spontaneous emission of 
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Figure 7: Left: Holon dispersion relation defined by Eq. (134). Right: Allowed 
values of q for adjacent odd N. Only the negative-energy holons are eigenstates 
of the Hamiltonian as positive-energy holons can lose energy by spontaneous 
emission of spinons. 



spinons. Let us assign to this holon the energy of the exact eigenstate onto 
which it projects and write 



pholon _ £ \1 

m - 2 \ N ) 



m (iV±l +m ) 



-E±l < m <0 
0< m <E=± 



(133) 



or 



E 



holo 



-E n 



(mod 2ir) 
(mod 2ir) 



(134) 



the crystal momentum q and spinon energy E q are defined as in Eqs. ( |43| ) and 
(^6|). This is shown in Fig. 7. The positive-energy part of the holon band is 
then unstable because it has negative curvature and thus does not satisfy the 
Landau criterion for spontaneous decay. A negative-energy holon, on the other 
hand, is forbidden from decaying by momentum conservation. 
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Figure 8: Contours used in Eqs. (136) and (137). 



A Fourier Sums 

Since the lattice sites z a are roots of unity we have 

N 

l[(z-z a )=z N -l . (135) 

a 

Then for < m < N we have 



and 



z a -l 



a=l 



z" 1 - 1 dz N f z m - x dz 



N r 

2rif c (z- l)(z N - 1) " 2m J c , (z - l)(z N - 1) 



1 ,m - 1, , ,m — 1, 2 

_iv + ( i + ( 2 dx 

2mf\,N, ,N. , N , , J i 2 

k ( ! ) +( 2 + ( 3 ) x* + ... } 



N + l . . 

— m , (I'JO) 



K-l m JV-1 m+ i 



^|z Q -l| 2 " ^-(z Q -l) 2 

a— 1 a— 1 ' 



n r z m dz n r z m dz 



2m Jc {z - l) 2 (z N - 1) 2m J c , (z - \) 2 {z N - 1) 
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1 



,AT. ,7V, ,AT. , J i 3 

( x ) +( 2 )x+( 3 )ar 2 + ... J 



7V 2 -1 m(N-l) m(m-l) 
12 2 + 2 



(137) 



3G 



B Problems 

1. Let r)i, rj N / 2 be the sites complementary to zi, z N / 2 - Show that 



N/2 N/2 N/2 N/2 

n - zfe ) 2 n z j = n ^ - ^ n n • 

j<k j j<k j 

2. Show that for any polynomial p(z) of degree less than N 



E 

P^a 



N P (0) - ( )p(z a ) 



3. Show that for a ^ /3 ^ 7 



[(4 • Sp), (S a x = 



[(S a -S J ),(S a xSp)]=i 
4. Show that the operator 



(S a ■ Sp)Sj — (Sj ■ S/3)S 



R a {z) 



(138) 



(139) 



(140) 



S a ■ (Sp x S 7 ) = 5 7 • (S a X Sp) = -S a ■ (S 7 X ^) (141) 



(142) 



(143) 



= E 



1 — z \ z a 



l + z\ Zp 



[i(S a xSp) + Sp] , (144) 



2 J z a - zp \ 2 J z a - zp 

where z is an arbitrary complex number, satisfies 

R a (z)\^>=0 , 
where \^f> is the Haldane-Shastry ground state, and 



(145) 
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JV-U, „ \ 1 
1 » 2 



S*- I £±£_ )5 • A 



(146) 



5. Show that the ground state wavefunction for a nonintcracting fcrmi sea 
on this lattice is 



N/2 N/2 

$(zi,...,z N /2,r)i,...,r) N/2 ) = Y[( z jVjy N/2 Y[( z j ~ z k)(Vj ~Vk) • (147) 



Then show that the Haldane-Shastry ground state is the Gutzwiller pro- 
jection of |$ >, i.e. that \^hs >— P \& >, where P is defined as in Eq. 
(107). Hint: If p denotes a permutation of N/2 things and sgn(p) is its 
sign, then 



(Af/2)! 

p 



N/2-1 
) X ■■• X z p (N/2) 



N/2 

Y[( Z J - z k) 

j<k 



(148) 



6. For lattice site z a different from 1 show that 



2 N (N ~ 1)/2 

— f—L\ 2 \ ^ pholon m 

j^tt' ^ m a 

m=-(7V+l)/2 



(N+l)/2 _ (N-l)/2 



(i - z a y 



l + z n 1 



I- z a 2 



z (N+l)/2 + z (N-l)/2 _ 2 



(149) 



where E%~ on is given by Eq. (|133|). Then use this result to show that 



<1p a \lp/3>= X^-X]^^ 1 ' ■- Z M)lpf3{zi, ...,%) 



2 <« a > f 1 ~ «^) (JV+1)/2 1 + (z*z p Y N -^ 2 



N { 1 - z* a z p 



(150) 



where \ip a > and > are localized spinon wavefunctions defined per Eq. 
4|). Hint: The fourier sum measures < ip a \Ca±l~LKYCpi\ipp >■ 
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7. The amplitudes at defined by Eq. (|5^) are the right eigenvectors of the 
non-hermitian matrix of Eq. (58). Show that the left eigenvectors of this 
matrix are 



2£[£ + m - 



1 

-n.1/2] £ (m-n+2k)b k 



(6o=0) , (151) 



k=l+l 



for £ < 0. The size of these "adjoint" coefficients may be judged from Fig. 
5, where we have plotted the function f m n{ z ) — Y]( biz 1 . Then show that 



N fe/2 
a j=Q 



(152) 



where \tp > is the Haldane-Shastry ground state \i\) m n > is the 2-spinon 
eigenstate defined by Eq. (|48|), and Cj are a set of coefficients. Obtain an 
expression these in terms of the bg. 



Let A/3 and Bp be any quantum operators and let 
a(z) 



b{z) = 



2N 



y N/2 _ -N/2 







Zp 



2N 



N 

£ 



z + zp 

Z - Z/3 



Show that 



(153) 



2iri 



(154) 



9. Show that the operator a(z) defined by Eq. (101) is hermitian whenever 
1*1 = 1- 

10. Let N be odd and let z\,...,Zm and rj\,...,r)M with M — (N — l)/2 be 
distinct lattice sites. Show that 



M 



N 2 -1 



j^k 1 IJ "" j^k 1 ■> n > j=l 1 31 

where h denotes the one site not occupied by Zj or rjj. 



12 



(155) 
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